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In order to obtain proper wave-function renormalization constants for unstable fermion and 
consist with Breit-Wigner formula in the resonant region, We have assumed an extension of the 
LSZ reduction formula for unstable fermion and adopted on-shell mass renormalization scheme in 
the framework of without field renormaization. The comparison of gauge dependence of physical 
amplitude between on-shell mass renormalization and complex-pole mass renormalization has been 
discussed. After obtaining the fermion wave-function renormalization constants, we extend them 
to two matrices in order to include the contributions of off-diagonal two-point diagrams at fermion 
external legs for the convenience of calculations of S-matrix elements. 
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I. INTRODUCTION 



We are interested in the LSZ reduction formula [1] in the case of unstable particle, for the purpose of obtaining a 
proper wave function renormalization constant (wrc.) for unstable fermion. In this process we have demanded the form 
of the unstable particle's propagator in the resonant region resembles the relativistic Breit-Wigner formula. Therefore 
the on-shell mass renormalization scheme is a suitable choice compared with the complex-pole mass renormalization 
scheme [2], as we will discuss in section 2. We also discard the field renormalization constants which have been proven 
at question in the case of having fermions mixing [3]. The layout of this article is as follows. First we discuss the 
problem arising from the introduction of field renormalization constants in the case of having fermions mixing. Next 
t-H ' we introduce the extended LSZ reduction formula for unstable fermion and obtain the proper fermion wrc. Next 
we discuss the difference between different renormalization schemes and the gauge dependence of physical amplitude. 
Next we change the fermion wrc. into two matrices for the purpose of including the corrections of fermion external 
, legs, thus simplify the calculations of S-matrix elements. Finally we give our conclusion. 

Along with the improvement of the precision of experiment, the radiative corrections need to be calculated more 
Oh! precisely, especially for the imaginary part of Feynman diagrams. For unstable particle the imaginary part of their 
self-energy functions should be considered in their wrc. [3]. Firstly let us introduce some notation in order to write 
them down. The bare fermion field is renormalized as follows 



m = Z^S>, ^ = ^Z^, (1) 

■ with the field renormalization constants and [3] 
H ' 

Zi = Z L ^ lL + Z R ^ lRl = Z L * lR + Z R * lL - (2) 

with 7£ and j R the left- and right-handed projection operator. Within one-loop accuracy we have Z^ = 1 + \5Z and 
2i = 1+ \5Z. Thus the inverse fermion propagator is as follows (the letters i,j denote the family indices) 

-iS^) = ^-m i )S ij +t ij {^), (3) 

with the one-loop renormalized self-energy 

± l3 (p) = ^(Sf ? (p 2 ) 7 L + ^(P 2 )1R) + m( P 2 )(m llL + m jlR )+ 



\5Z i:j (^ - rrij) + \(j/>- mi)SZij - 8mA 



(4) 



If we introduce the field renormalization conditions that the fermion two-point functions are diagonal when the 
external-line momcntums are on their mass shell, 



S lJ 1 (p)Uj(p)\ p 2 =m 2 = 
^(p')S~ j l {p')\ pl 2 =m 2 = 



i*3- (5) 
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with Ui and Uj the Dirac spinors, we will encounter a bad thing that there is no solution for the above equations if 
we keep the " pseudo-hermiticity" relationship between zi and Zi [3] 

Z3=7°Z'V (6) 

This is because of the existence of the imaginary part of the fermion two-point functions, which makes E^- , E^ and E?- 
un-Hermitian [4]. As we know the relationship of Eq.(6) comes from the hermiticity of Lagrangian, which guarantees 
the S-matrix is unitary, thus cannot be broken. 

A method to solve this problem is to introduce two sets of field renormalization constants, one set for the external- 
line particles, the other set for the inner-line particles in a Feynman diagram [3]. But it brings some complexity and 
cannot explain why we couldn't find a consistent field renormalization constant for all particles. 

II. WAVE-FUNCTION RENORMALIZATION CONSTANTS 

Another method to solve the above problem is to abandon field renormalization, since particle fields don't belong to 
physical parameters thus don't need to be renormalized. One effect of this method is that the calculations of Feynman 
diagrams becomes simpler than before. The other effect will be seen in section 3. 

Since there aren't field renormalization constants, we need to find a way to determine the fermion wrc. We demand 
the propagator of unstable particle has the following form when the momentum is on mass shell 



p 2 — to 2 + imT(p 2 ) 



(7) 



where mT{p 2 ) is a real number which represents the imaginary part of the self-energy function times the particle's 
wrc. This form resembles the relativistic Breit-wigner formula, thus will be a reasonable assumptions. Under 
this assumption only the on-shcll mass renormalization scheme is suitable, the other scheme, complex-pole mass 
renormalization scheme, cannot satisfy Eq.(7). We can see this point in the boson case. 

The mass definition in complex-pole scheme is based on the complex- valued position of the propagator's pole [2]: 

s = M 2 + A(s) . (8) 

with A(p 2 ) the self-energy function of boson. Writing s — to 2 — im 2 r 2 , the mass and width of the unstable particle 
may be identified with to 2 and T 2 , respectively. Given to 2 and r 2 , there are two other definitions: 

mi = ^/mi + r 2 ,r 1 = ^r 2 , 

s = (m 3 - i^) 2 . 
The form of the propagator in resonant region is 

1 1 _ (l-^m 2 ))- 1 

p 2 - M 2 - A(p 2 ) ~ m 2 - M 2 - A(m 2 ) + (p 2 - to 2 )(1 - A'(m 2 )) ~ p 2 - m 2 + B 

with B = (to 2 — Mq — A(to 2 ))/(1 — A'{m 2 )). Through simple calculations we find all of these three definitions of 
mass cannot make B a pure imaginary number. So the complex-pole mass renormalization scheme isn't suitable for 
the assumption of Eq.(7). 

Because the unstable fermion's propagator near resonant region must contain Dirac matrices, Eq.(7) isn't enough. 
We further assume it has the following form: 

i p 2 -> to 2 iZf {j) + rrn + ix)Zf ^ 



j)-m t + ± u {i)) p 2 -ml + imiTdp 2 )' 

i _ i 

where Zf and Z? are the fermion wrc, a; is a quantity which has nothing to do with our present calculations. The 
LSZ reduction formula which relates the Green functions to S-matrix elements may have such form: 
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-2 _^ „_,2 



/ dxfdxje ipiXl e ip2X2 < fl\Tip(x 1 )ip(x 2 ) ■ ■ ■ \Sl > 

(10) 



Pi ~> W ? Zfjfi+mi+ixi) 7/ \ . . . I CI . . . 1 > (]>2+^+^)^ 



here we have assumed the Heisenberg field i[>(xi) can be transformed into the out-state <ip(pi)\ with an additional 
i 

coefficient pfi^+imir"^) ' an ^ tnc Heisenberg field ip(x2) can be transformed into the in-state \ij){pz)> with an 
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additional coefficient p ^^"^rJ^) • According to Eq.(9), we can sum up the corrections to each external legs in a 
Green function and analyse the Green function into the amputated part times the corresponding propagators. That 
is to say 

/ dxfdxje ipiXl e ip2X2 < Q\T-4){x 1 )ip{x 2 ) ■ ■ ■ |fi > 

where the subscript amp refers to the amputated Green function. Thus the LSZ reduction formula with two fermions 
at external legs will be: 



<p 1 ...p n \S\k 1 ...k m >= Z^z—u{ Pl )Zf (Ampj Z;u(fci). (12) 

k\ k m 

_ i i 

with Zf and ZJ the fermion wrc. and Z the boson wrc. The fermion spinors u(pi) and u(k\) have been added to 
the right hand side of the above equation, as are demanded. 

In actual calculations, the fermion propagator needs to be transformed into the form of Eq.(9). To start with, we 
introduce the renormalized fermion self-energy functions as follows 

s« = i>(^{p 2 YiL + sgo^hfl) + e£ Vhi + s« Vh* - . (i3) 

Although the self-energy functions contain only the contributions of one-particle-irreducible (1PI) diagrams in usual 
sense, it isn't denied in principle to include the contributions of reducible diagrams in the self-energy functions. In 
the present scheme the non-diagonal two-point functions aren't zero when external momentum is on shell, since there 
aren't the conditions of Eqs.(5). Thus the diagonal self-energy functions will contain the contributions of the reducible 
diagrams: 




k,l,- ■ -,171 

k, I, ■ ■ -,m 7^ i 



(14) 



where the disc is the 1PI diagram. The condition k,l,- ■ -,m ^ i avoids the reduplication in actual calculations. Now 
we write the inverse fermion propagator as follows 



with 



-iSii 1 ^) =p'-m l + tu(p) = j>{a- lL + bj R ) + cj L + dj R 

a= 1 + Sf 4 (p 2 ), b= 1 + S«(p 2 ), 

and therefore after some algebra [3] 



c = Sf/ L (p 2 ) - mi -5 mil d = ^fi R (p 2 ) -mi- 5m t . 



i>{a-y L + bj R ) - dj L - cj R 
p 2 ab — cd 



(15) 
(16) 

(17) 



That's to say 
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i(rm + 5m l+ jr iL (l + ££) + fr R (l + Eg) - 7L £^ L - 7K £?' fl ) 



(18) 



ji - tth + S«(|J) ~ ' P 2 (l + ££)(1 + Ej?) - (m 4 + «5m 4 - Ef^K + 5m { 
Adopting the on-shell mass renormalization scheme 1 , the mass renormalization condition is 

Re[m 2 (l + ££(m?))(l + ££(m?)) - K + Sm t - £*< L (m?))K + 5m t - Ef^ro?))] = (19) 
Expanding the denominator of Eq.(18) at p 2 = mf according to Eq.(9), we obtain 



j P Z ^ m j *K + gmj +pj L (l + g|(m|)) + ^(1 + Eff(m 2 )) - 7Lg|V|) ~ 7gg|Vf)) f20) 

t-mi + tuU) {p 2 -m 2 i +im i T i {p 2 ))A 



with 



A = Re[l + E£ + S|? + ££s£ + m|^(Ef 4 + E£ + ££e£)+ 
^gi^ + ^-Ef^ + ^K+Jm.-Ef/ 1 ')] 

m^ip 2 ) = l/mb 2 (Ef(/) + E«(p 2 ) + E^ 2 )E«( P 2 )) + (m, + fm,)(E* V) + ^ V)) - V)E& V)] 



(22) 



Noted we have discarded the regular terms in Eq.(20). 
From Eq.(9) and (20) we obtain 

Z-\Zi \ = (1 + Si?(m?))/A, zf*Zf* = (1 + ££(mf))M, 

zf*zf* - (m» + frm - ))/(^K + ix)) , Zf h 2f ' - (m l + 5™ 4 -S^ L (m 2 ))/(A(m 4 + za;)). 



(23) 



In order to solve these equations, we firstly obtain a equation about x 

( mi + ix) 2 (l + ££(m?))(l + £g(ro?)) = (m, + im, - E^ L (m 2 ))(m, + ,5m, - Ef^(m 2 )) ■ (24) 
Replacing + z.t in Eqs.(23) by this relationship, we obtain 

= (m,+3m 8 -S| fl K 2 ))(l+SgK 2 )) 
(m i +«m i -E^( m 2))(l+E l f i ( m 2)) ' 

(m 1 +5m l -S^ K (mf))(l+Ef i (m2)) ' 

Because we haven't introduced the fermion field renormalization constants, which contain Dirac matrices thus may 
change the Lorentz structure of the fermion's self-energy functions, the self-energy functions in our scheme will keep 
the Lorentz invariant structure: 

%j (f) = m% (P 2 )1L + Eg (p 2 ) lR ) + Eg (p 2 ) (m ,, 7 L + moj7fl) - Wtf • (26) 
where mo,» and m 0j the bare fermion masses. It is obvious that 

= Ef/ L = m.E? (27) 

Therefore Eqs.(25) can be simplified as follows 



1 There are many articles discussing the gauge dependence of the on-shell mass renormalization scheme [5], we will discuss this 
problem in section 4 
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Z^Z^ = (l + ££(m?))M Zf^zf^ = (1 + E£(m?))/A, (28) 

Zf/Zf =(l + S>, 2 ))/(l + S^K)),f/^ = (1 + Eg(m?))/(1 + E£(m?)). 

Because the four equations aren't independent, we need to find an additional condition. From the strong symmetry 
between Z and Z manifested in the above equations, a reasonable assumption is to set Zf = Zf , Zf = Zf, which 
is also the case of stable fcrmion. Therefore the final results are 



Zf = Zf = (1 + Eg(m?))/A, 

zf = zf = (i + Ef i K))M. 



(29) 



and 

Srrii 
mi 



= ^i?e[(l + E£(m?))(l + E«(mf ))] + Jm[E«(m?)P + Re^rnf)] 1 (30) 



It is very easy to give the explicit one-loop level results. The wrc, mass counterterm and decay rate are listed 
below: 

Zf = Zf = l- ifep&ro?) + mf ^(Ef 4 (p 2 ) + E«(p 2 ) + 2Ef 4 ( P 2 ))| p2=m2 ] + i/m[E£(m 2 )] , 

= Zf = 1 - ite[Ej?(m?) + m 2 ^(E^ 2 ) + E«(p 2 ) + 2£f>(p 2 ))| p2= „; 2 ] + J /m[£f 4 (m 2 )] . ^ 

Jm, = ^i? e [E£(m 2 ) + E«(m 2 ) + 2£i(m 2 )] (32) 

™j, - m 2 /m[£ 2 L 4 (m 2 ) + E«(m 2 ) + 2£f i (m 2 )] (33) 

Another meaningful result concerning x (which appears in Eq.(9)) is: 

x = -Ti/2 (34) 

Until now we haven't discussed the Optical Theorem. It is a very complex problem, so we want to leave it for the 
future work. Here we only point out that Eq.(33) is really the fermion decay rate at one- loop level [3]. 

III. COMPARISON BETWEEN DIFFERENT RENORMALIZATION SCHEMES 

Adopting which mass renormalization scheme and whether introducing field renormalization constants determine 
the difference of the renormalization schemes. Here we will list three other schemes and compare their results. 

Firstly we consider the renormalization scheme which adopts the complex-pole mass renormalization scheme and 
doesn't introduce field renormalization constants. In this case Eq.(9) isn't satisfied, but we can assume Tj is a 
complex number and independent of the propagator's momentum. Thus all of the results associated with fermion 
wrc. in previous section keep unchanged, except for removing the operator Re from A, because in this case it is no 
sense to treat the real part and the imaginary part of the denominator of the fermion propagator differently. Therefore 
we obtain 



Zf = Zf = (l + EiKm?))/^ 
Zf = Zf = (l + E£K))Mi 



(35) 



with 



A 1 = 1 + ££ + £* + Ef 4 £« + m\—^l + £« + e£e«) + -^-{rm + 5m t - £*•*) + -^-(m, + Srm - 

(36) 

On the other hand, Sm,i and Ti arc different from the previous results very much. But at one-loop level they are the 
same as Eq.(32) and (33). 

Next we study the renormalization scheme which adopts the complex-pole mass renormalization scheme and has 
introduced field renormalization constants. We still assume Ti is a complex number and independent of the propaga- 
tor's momentum. Because the field renormalization constants have been introduced, the fcrmion self-energy functions 
will contain their contributions. We can write down the fermion self-energy functions in such form: 
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= £&*7 L + Ei?M + S^Sl +4' • (37) 



-iS,.R 



where x denotes the field renormalization constants. This form can be easily understood if we regard it as an 
application of the LSZ reduction formula Eq.(12) at two-point Green functions. The new self-energy function £' is 
called "pure" self-energy function by us. It is easy to obtain from the above equation that 

y>L _ y L k y L k^piL y-R _ y^k 7^5y/fl 

'-'ii ~ ^ii ^ii ^ii i ^ii ~ n ii ^ii '-'ii ' /oo\ 

yS,L y^2 2 -ytS .L ir^S,R, y^2 y^2 y'S.R 

^ii — ^ii ^ii ^ii ! ^ii — ^li a ii ^ii 

We can write down the inverse fermion propagator in the form of Eq.(15) with 

« = ZpZp + *h(p 2 ) , b = zfzf + S«(p 2 ) , 

c = Ef.'Vj-K + ^O^^.d = Sfi'Vj-K + ^mi)^^'. 

According to Eq.(17), the condition that in the limit p 2 — > mf the chiral structure in the numerator must be canceled 
out leads to 



(40) 



£*• V? ) - K + SmJZpzf:* = V?) - (m 4 + tfmO^Z?* . 
On the other hand the unit residue renormalization condition amounts to requiring 

a = ab + ™ 2 a-^b + rajb^a - c^d - d^c (41) 

Using Eq.(39) and Eqs.(38), we have 

1 _ y L i y L i ri i v'i _i_ ™2/ a V '£ _i_ j±gil 9 V'^ j_ m '+ l5m '~ s »f'" n 9 vVS,L roijjmi-S^ g V /S,fli 

1 - ^ii ^li i 1 + ^ii + m i Kdp^^ii + i+Ejf dp* ™ ' H SF^ii H " J ^ ' 

Like the case in section 2 we need an additional condition to determine the field renormalization constants. If we 
adopt the same assumption as the ones in the previous section, Zfi = Zfl , Z^ = Z^, we will obtain 



Zh = Z& = (l + E{f(m?))M 2) 
Z% = Zn = (1 + I#(m?))/A 2 . 



(43) 



with 



ii ii ' i dp' 2 V ii ii ' ii ii 



' •■■ . , (44) 



Ai is the same as A\, except for replacing the self-energy functions with "pure" self-energy functions. But at one-loop 
level they are the same, so the one- loop results of Eqs.(43) are [3] 

Zt =Zl = \- ££(m?) - m 2 ^(£f 4 (p 2 ) + S«(p 2 ) + 2£f 4 (p 2 ))| p2=m2 , 

Zf = Zf = \- Ei|(m 2 ) - m?^(Ef i (p 2 ) + E«(p 2 ) + 2^))\ v 2^2 . ^> 

Finally we consider the renormalization scheme which adopts the on-shell mass renormalization scheme and has 
introduced field renormalization constants. In this scheme, the self-energy functions can still be written down in the 
form of Eqs.(37) and the renormalization conditions of Eqs.(40) keep unchanged. The only variation comes from the 
unit residue condition. In on-shell scheme we should treat the real part and the imaginary part of the denominator 
of the fermion propagator in different ways. Thus the unit residue condition becomes 

d d d d 

a = Re[ab + mja-^b + m 2 b^a - c^d - d^c] (46) 

after some algebra we have 
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mj+faii-g^ g y /S,L m t +Sm t -T,[f L g wS.iiii 
(l+S^)(l+S^f)9F « + (l+E£)(l+£<f) JJ 

If we adopt the consistent assumption: Zf — Zf- , Z-f — Zf, we will have 

Zh = Zk = 1/(4,(1 + E£(m?))), 
ZH = zn = 1/(4(1 + £<?(m?))). 



(47) 



(48) 



with 

2 

3 +H 'h + md P i ^ + i + v£d P i « + (i + i#)(i + £<fW « + (i + E^)(i + s;f)ap2 



</ inj u v// - + <5to z - Ej^ 9 /S|L to, + <5to z - gjj ' L (9 /S)fl 



(49) 



(50) 



which arc different from all of the previous results. At one- loop level, Eqs.(48) becomes 

Zt = Z}- = 1- ££(m 2 ) - m^Re^p 2 ) + S^(p 2 ) + 2£g(p 2 )] pa=in j , 
Zf = Z? = 1- E£(m?) - mf^epSb 2 ) + S||(p 2 ) + 2£f 4 (p 2 )] p2=m , . 

which are different from Eqs.(31) and (45). 

Through these comparison we find that whether to introduce field rcnormalization constants or not may change the 
fermion wrc. (comparing Eqs.(50) with Eqs.(31)). On the other hand, introducing field renormalization constants will 
enhance the complexity of calculating S-matrix elements. As it's very known, "bare" field isn't a physical quantity 
thus doesn't need to be renormalized. So the best way to renormalize a quantum field theory is to discard field 
renormalization. It not only simplifies the calculations of S-matrix elements, but also avoids the possible deviation 
arising from the introduction of field renormalization constants. 



IV. GAUGE DEPENDENCE OF RENORMALIZATION SCHEMES 



In this section we will discuss the gauge dependence of the rcnormalization schemes we have mentioned. For 
convenience we only discuss two typical renormalization schemes. The first scheme is the scheme we have adopted, 
the second scheme is the scheme which adopts the complex-pole scheme and contains field renormalization constants. 
At one-loop level, the results of Srrii and Ti are the same in these two schemes. The difference comes from the imaginary 
part of the fermion's wrc, if you compare Eqs.(31) with Eq.(45). So we only need to consider the imaginary part of 
the fermion self-energy functions in the following calculations. 

We have selected two physical processes to test the gauge dependence of the physical amplitude in these two 
renormalization schemes. One process is W gauge boson decays into a lepton and an anti-neutrino. For our purpose 
it is sufficient to only consider the dependence of W boson gauge parameter £w Thus there are only two Feynman 
diagrams which contribute to the imaginary part of the lepton self-energy functions, as shown in Fig.l 



G W 




V; V; 

FIG. 1. lepton's self-energy diagrams which have contributions to the ^-dependent imaginary part of lepton's self-energy 
functions 



Using the cutting rules [6], we can calculate the imaginary part of the lepton's wrc. In the first scheme we have: 
(restrained > 0) 

'™t Z < 1 " 64»MS,» W , < 51 » 

with mi y i the lepton's mass, Mw the W gauge boson mass, sw the sine of the weak mixing angle 6w and 9 the 
Heaviside function. 
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FIG. 2. irreducible one-loop diagrams for W — ► ejf, which have contributions to the £w-dependent imaginary part of this 
amplitude 



In Fig. 2 we show the irreducible diagrams that contribute to the ^-dependent imaginary part of the one-loop 
W~ — > CiVi amplitude. Using the cutting rules we obtain 



with 



and the form factors 



F L 



Jm[Ti] - A Lj F 1 L +B L G L +^ L /m[Zf]e/(2V2 SH ,) 
-^l = u{p 2 )ilLv{vi) ,B L = u{p 2 )^~i L v{pi) ■ 



48v / 2u's H / 
e 3 (fa-l) 3 "i|,iM ff 8[l-^] 
32v / 27rs w (M^,-m^ i ) 



64v / 27rM^s^mf ! , 



(52) 
(53) 

(54) 



Thus we have 



D e J (g 1 y-l) 3 mi, i M w 9[l-gw] 
L 32V2irs w (M^ v -mf_ i ) 



(55) 



We find that the imaginary part of the one-loop W — ► e^i amplitude is gauge dependent in our scheme. On the 
other hand, the second scheme gives 



Im[Zt] = 



L _e\ml i -M^ w fe[ml i -M^ w ] 



and the imaginary part of the one- loop W — ► amplitude (according to the above results) is 



Im[T{] 



-A, 



° 3 (&~3& + 2)0[l-fr y ] , _ e 3 (^-l) 3 m M M^[l-^] 



48V2irs w 



+ B 



32V2irs w (M^ -mf t ) 



(56) 



(57) 



also a gauge dependent quantity. At one-loop level Eq.(55) and (57) don't contribute to the modulus square of the 
physical amplitude. But at two-loop level they will have the contributions: Im[Ti^ Ira\Ti] or Im[T{]i Im[T[\. At the 
limit mis — * they have the same contributions as follows 



JmpTilt/mpTiU^o = Im[T^ Im[T[] mi .^ 



2304tt 2 s^ 



(58) 
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Of course there will be no such gauge-dependent terms if we adopt the unitary gauge. We hope that the contributions 
from the two-loop diagrams of the process W~ — -> e^P, may cancel these gauge-dependent terms. 

Another process we have studied is Z gauge boson decays into a pair down-type quarks. The calculations are more 
complex than the previous ones. The imaginary parts of the down-type quark's wrc. are the same in these two 
schemes in the limit m^i — * (m^, is the down- type quark's mass). They are as follows 



Im[M(Z - di di)] md ^ = A L [- 384 \ » e 3 (l - Ac 2 w t w f/ 2 6[l - Ac 2 w t w ]- 



192 i aw <?((Zw - l) 2 c 4 w - 2(£ W - 5)c 2 w + l)y/(£ w - l) 2 c* w - 2(t w + l)c 2 w + 10[l/c w - - 1]] 
The contribution of this quantity to the modulus square of the Z — > didi amplitude at two- loop level is 
Im[M(Z -+ didi)]1lm[M{Z - <^)] m<M ^o = ~ 7372 &X» M^c^w - \fQ\\ - 4c 2 w ^ w }+ 

WW 

? J cH gi e 6 M^((g w - l) 2 c A w - 2(£ w - 5)c 2 w + l)9[l/c w - - 1] 



18432,, ^ w ° w 

(switiw - l) 4 c^ - 4(fry - 2)(fry - l) 2 c 6 w + (6^ - 20£w - 18)4, - 4 (6v - 2)<& + 1)- 
(1 - 4c^) 3/ V(£w - \) 2 c* w - 2(t w + l)c 2 w + 1) 



(59) 



(60) 



There are also gauge-dependent terms. But we don't know how they can be cancelled in the full modulus square of 
the Z — ► didi amplitude at two-loop level. This problem will be left for our next work. 

Through these discussions we can draw a conclusion that based on the present knowledge we cannot judge which 
scheme is the better. 



V. INTRODUCTION OF OFF-DIAGONAL WAVE-FUNCTION RENORMALIZATION CONSTANTS 

In order to simplify the calculations of the corrections of fermion external legs in a Feynman diagram we can change 
the fermion wrc. to two matrices, which will contain the contributions of the off-diagonal two-point diagrams at 
fermion external legs. In order to simplify this definition we need to define two new self-energy functions as follows: 



(61) 



k,l,- ■ -,m 7^ i 

and the corresponding renormalized self-energy functions (after adding proper counterterms to the above diagrams): 





(62) 



We will see in a second that the new self-energy functions have no other meanings except for bringing about great 
convenience in actual calculations. 

_ i i 

Now we introduce the fermion wrc. matrices Zf. and Zf- as follows: 



<Pi...\S\kj... >= u(m l )Z l 2 MijZ?u(mj) 



«(mo4^7 z *i u K-). 



(63) 



with Si^j = 1 when i ^ j and M amp the amplitude of completely amputated diagrams (without off-diagonal self- 
energy diagrams at fermion external legs). Noted we have plainly used the definitions of Eqs.(62) in the second line 
of Eqs.(63). The last line of Eqs.(63) gives the definitions of fermion wrc. matrices: 



Zluimj) = {5ij + J^iXi^Zjuim^ 
i 



u(mi)Zl = u(m i )Z?(a y +*EyG} i ) ? ^)| p?= 



(64) 



Pt= m i ' 
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after some algebra we have 



Zf = Zfhij + j^fam^f *E£.(m?) + m*Z? *£«(m?) + m^EgV?) + ro^f *Eg fl (m?)] , 
= Zfh l0 + ^[m imj Z^{rn*) +m *zf^^ . 



At one-loop level these results agree with Eqs.(3.3) and Eqs.(3.4) in Ref. [3]. 

Using these wrc. matrices we can calculate S-matrix elements without considering any of the corrections of fermion 
external legs in a Feynman diagram. Noted such matrices are also suitable for gauge bosons which have mixing 
between each other. 



VI. CONCLUSION 



In this paper we have attempted to extend LSZ reduction formula for unstable fermion. We adopted the assumption 
of Eq.(7), which resembles the relativistic Breit-Wigner formula thus can describe the behavior of unstable particle's 
propagator in the resonant region. Under this assumption only the on-shcll mass renormalization scheme is suitable 
(see section 2). In consideration of the difficulty of introducing fermion field renormalization constants [3] we have 
abandoned field renormalization. We assume the unstable fermion's propagator has the form of Eq.(9) in the limit 
p 2 — ► mf. Thus we can get the LSZ reduction formula Eq.(12) for unstable fermion and obtain the proper fermion 
wrc. Eqs.(29). 

Next we change our assumption and compare the difference between different renormalization schemes. There are 
four basic renormalization schemes, corresponding to whether or not to introduce field renormalization constants 
and to adopt which mass renormalization scheme, on-shcll scheme or complex-pole scheme [2]. Through comparison 
we find that whether to introduce field renormalization constants or not will affect the values of fermion wrc. and 
demonstrate its difference in the imaginary part of physical amplitude, and so does the choice of mass renormalization 
scheme (see section 3 and 4). This discrepancy is gauge dependent and can manifest itself at the modulus square of 
physical amplitude at two-loop level. But we haven't found strong evidence to judge which scheme is better that the 
others. 

Next we change the fermion wrc. to two matrices, in order to contain the corrections of fermion external legs in a 
Feynman diagram. The new fermion wrc. matrices, according to Eqs.(63), will simplify the calculations of physical 
amplitude, because we don't need to calculate the corrections of off-diagonal self-energy diagrams at fermion external 
legs. We note that the concept of wrc. Matrix can be extended to bosons provided they have mixing in themselves. 

In conclusion, we recommend the renormalization scheme which adopts the on-shell mass renormalization scheme 
and abandons field renormalization, because this renormalization scheme is very simple and clear. We will study this 
point further in the future. 
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